WHO elements - A new category of selfish genetic elements at the borderline between homing elements and transposable elements.

Homing genetic elements are a form of selfish DNA that inserts into a specific target site in the genome and spreads through the population by a process of biased inheritance. Two well-known types of homing element, called inteins and homing introns, were discovered decades ago. In this review we describe WHO elements, a newly discovered type of homing element that constitutes a distinct third category but is rare, having been found only in a few yeast species so far. WHO elements are inferred to spread using the same molecular homing mechanism as inteins and introns: they encode a site-specific endonuclease that cleaves the genome at the target site, making a DNA break that is subsequently repaired by copying the element. For most WHO elements, the target site is in the glycolytic gene FBA1. WHO elements differ from inteins and homing introns in two fundamental ways: they do not interrupt their host gene (FBA1), and they occur in clusters. The clusters were formed by successive integrations of different WHO elements into the FBA1 locus, the result of an 'arms race' between the endonuclease and its target site. We also describe one family of WHO elements (WHO10) that is no longer specifically associated with the FBA1 locus and instead appears to have become transposable, inserting at random genomic sites in Torulaspora globosa with up to 26 copies per strain. The WHO family of elements is therefore at the borderline between homing genetic elements and transposable elements.

The Transcriptional Regulator Prdm1 Is Essential for the Early Development of the Sensory Whisker Follicle and Is Linked to the Beta-Catenin First Dermal Signal.

Prdm1 mutant mice are one of the rare mutant strains that do not develop whisker hair follicles while still displaying a pelage. Here, we show that Prdm1 is expressed at the earliest stage of whisker development in clusters of mesenchymal cells before placode formation. Its conditional knockout in the murine soma leads to the loss of expression of Bmp2, Shh, Bmp4, Krt17, Edar, and Gli1, though leaving the ß-catenin-driven first dermal signal intact. Furthermore, we show that Prdm1 expressing cells not only act as a signaling center but also as a multipotent progenitor population contributing to the several lineages of the adult whisker. We confirm by genetic ablation experiments that the absence of macro vibrissae reverberates on the organization of nerve wiring in the mystacial pads and leads to the reorganization of the barrel cortex. We demonstrate that Lef1 acts upstream of Prdm1 and identify a primate-specific deletion of a Lef1 enhancer named Leaf. This loss may have been significant in the evolutionary process, leading to the progressive defunctionalization and disappearance of vibrissae in primates.

Exploiting epigenetic dependencies in ovarian cancer therapy.

Ovarian cancer therapy has remained fundamentally unchanged for 50 years, with surgery and chemotherapy still the frontline treatments. Typically asymptomatic until advanced stages, ovarian cancer is known as "the silent killer." Consequently, it has one of the worst 5-year survival rates, as low as 30%. The most frequent driver mutations are found in well-defined tumor suppressors, such as p53 and BRCA1/2. In recent years, it has become clear that, like the majority of other cancers, many epigenetic regulators are altered in ovarian cancer, including EZH2, SMARCA2/4 and ARID1A. Disruption of epigenetic regulators often leads to loss of transcriptional control, aberrant cell fate trajectories and disruption of senescence, apoptotic and proliferation pathways. These mitotically inherited epigenetic alterations are particularly promising targets for therapy as they are largely reversible. Consequently, many drugs targeting chromatin modifiers and other epigenetic regulators are at various stages of clinical trials for other cancers. Understanding the mechanisms by which ovarian cancer-specific epigenetic processes are disrupted in patients can allow for informed targeting of epigenetic pathways tailored for each patient. In recent years, there have been groundbreaking new advances in disease modeling through ovarian cancer organoids; these models, alongside single-cell transcriptomic and epigenomic technologies, allow the elucidation of the epigenetic pathways deregulated in ovarian cancer. As a result, ovarian cancer therapy may finally be ready to advance to next-generation treatments. Here, we review the major developments in ovarian cancer, including genetics, model systems and technologies available for their study and the implications of applying epigenetic therapies to ovarian cancer.

Giant GAL gene clusters for the melibiose-galactose pathway in Torulaspora.

In many yeast species, the three genes at the centre of the galactose catabolism pathway, GAL1, GAL10 and GAL7, are neighbours in the genome and form a metabolic gene cluster. We report here that some yeast strains in the genus Torulaspora have much larger GAL clusters that include genes for melibiase (MEL1), galactose permease (GAL2), glucose transporter (HGT1), phosphoglucomutase (PGM1) and the transcription factor GAL4, in addition to GAL1, GAL10, and GAL7. Together, these eight genes encode almost all the steps in the pathway for catabolism of extracellular melibiose (a disaccharide of galactose and glucose). We show that a progenitor 5-gene cluster containing GAL 7-1-10-4-2 was likely present in the common ancestor of Torulaspora and Zygotorulaspora. It added PGM1 and MEL1 in the ancestor of most Torulaspora species. It underwent further expansion in the T. pretoriensis clade, involving the fusion of three progenitor clusters in tandem and the gain of HGT1. These giant GAL clusters are highly polymorphic in structure, and subject to horizontal transfers, pseudogenization and gene losses. We identify recent horizontal transfers of complete GAL clusters from T. franciscae into one strain of T. delbrueckii, and from a relative of T. maleeae into one strain of T. globosa. The variability and dynamic evolution of GAL clusters in Torulaspora indicates that there is strong natural selection on the GAL pathway in this genus.

Polymorphic centromere locations in the pathogenic yeast Candida parapsilosis.

Centromeres pose an evolutionary paradox: strongly conserved in function but rapidly changing in sequence and structure. However, in the absence of damage, centromere locations are usually conserved within a species. We report here that isolates of the pathogenic yeast species Candida parapsilosis show within-species polymorphism for the location of centromeres on two of its eight chromosomes. Its old centromeres have an inverted-repeat (IR) structure, whereas its new centromeres have no obvious structural features but are located within 30 kb of the old site. Centromeres can therefore move naturally from one chromosomal site to another, apparently spontaneously and in the absence of any significant changes in DNA sequence. Our observations are consistent with a model in which all centromeres are genetically determined, such as by the presence of short or long IRs or by the ability to form cruciforms. We also find that centromeres have been hotspots for genomic rearrangements in the C. parapsilosis clade.

The yeast mating-type switching endonuclease HO is a domesticated member of an unorthodox homing genetic element family.

The mating-type switching endonuclease HO plays a central role in the natural life cycle of Saccharomyces cerevisiae, but its evolutionary origin is unknown. HO is a recent addition to yeast genomes, present in only a few genera close to Saccharomyces. Here we show that HO is structurally and phylogenetically related to a family of unorthodox homing genetic elements found in Torulaspora and Lachancea yeasts. These WHO elements home into the aldolase gene FBA1, replacing its 3' end each time they integrate. They resemble inteins but they operate by a different mechanism that does not require protein splicing. We show that a WHO protein cleaves Torulaspora delbrueckii FBA1 efficiently and in an allele-specific manner, leading to DNA repair by gene conversion or NHEJ. The DNA rearrangement steps during WHO element homing are very similar to those during mating-type switching, and indicate that HO is a domesticated WHO-like element.

In the same way as a sperm from a male and an egg from a female join together to form an embryo in most animals, yeast cells have two sexes that coordinate how they reproduce. These are called "mating types" and, rather than male or female, an individual yeast cell can either be mating type "a" or "alpha". Every yeast cell contains the genes for both mating types, and each cell's mating type is determined by which of those genes it has active. Only one mating type gene can be 'on' at a time, but some yeast species can swap mating type on demand by switching the corresponding genes 'on' or 'off'. This switch is unusual. Rather than simply activate one of the genes it already has, the yeast cell keeps an inactive version of each mating type gene tucked away, makes a copy of the gene it wants to be active and pastes that copy into a different location in its genome. To do all of this yeast need another gene called HO. This gene codes for an enzyme that cuts the DNA at the location of the active mating type gene. This makes an opening that allows the cell to replace the 'a' gene with the 'alpha' gene, or vice versa. This system allows yeast cells to continue mating even if all the cells in a colony start off as the same mating type. But, cutting into the DNA is risky, and can damage the health of the cell. So, why did yeast cells evolve a system that could cause them harm? To find out where the HO gene came from, Coughlan et al. searched through all the available genomes from yeast species for other genes with similar sequences and identified a cluster which they nicknamed "weird HO" genes, or WHO genes for short. Testing these genes revealed that they also code for enzymes that make cuts in the yeast genome, but the way the cell repairs the cuts is different. The WHO genes are jumping genes. When the enzyme encoded by a WHO gene makes a cut in the genome, the yeast cell copies the gene into the gap, allowing the gene to 'jump' from one part of the genome to another. It is possible that this was the starting point for the evolution of the HO gene. Changes to a WHO gene could have allowed it to cut into the mating type region of the yeast genome, giving the yeast an opportunity to 'domesticate' it. Over time, the yeast cell stopped the WHO gene from jumping into the gap and started using the cut to change its mating type. Understanding how cells adapt genes for different purposes is a key question in evolutionary biology. There are many other examples of domesticated jumping genes in other organisms, including in the human immune system. Understanding the evolution of HO not only sheds light on how yeast mating type switching evolved, but on how other species might harness and adapt their genes.

Genomic diversity and meiotic recombination among isolates of the biotech yeast Komagataella phaffii (Pichia pastoris).

BACKGROUND: Komagataella phaffii is a yeast widely used in the pharmaceutical and biotechnology industries, and is one of the two species that were previously called Pichia pastoris. However, almost all laboratory work on K. phaffii has utilized strains derived from a single natural isolate, CBS7435. There is little information about the sequence diversity of K. phaffii or the genetic properties of this species. RESULTS: We sequenced the genomes of all the known isolates of K. phaffii. We made a genetic cross between derivatives of two isolates that differ at 44,000 single nucleotide polymorphism sites, and used this cross to analyze the rate and landscape of meiotic recombination. We conducted tetrad analysis by making use of the property that K. phaffii haploids do not mate in rich media, which enabled us to isolate and sequence the four types of haploid cell that are present in the colony that forms when a tetra-type ascus germinates. CONCLUSIONS: We found that only four distinct natural isolates of K. phaffii exist in public yeast culture collections. The meiotic recombination rate in K. phaffii is approximately 3.5 times lower than in Saccharomyces cerevisiae, with an average of 25 crossovers per meiosis. Recombination is suppressed, and genetic diversity among natural isolates is low, in a region around centromeres that is much larger than the centromeres themselves. Our work lays a foundation for future quantitative trait locus analysis in K. phaffii.

Coverage-Versus-Length Plots, a Simple Quality Control Step for de Novo Yeast Genome Sequence Assemblies.

Illumina sequencing has revolutionized yeast genomics, with prices for commercial draft genome sequencing now below $200. The popular SPAdes assembler makes it simple to generate a de novo genome assembly for any yeast species. However, whereas making genome assemblies has become routine, understanding what they contain is still challenging. Here, we show how graphing the information that SPAdes provides about the length and coverage of each scaffold can be used to investigate the nature of an assembly, and to diagnose possible problems. Scaffolds derived from mitochondrial DNA, ribosomal DNA, and yeast plasmids can be identified by their high coverage. Contaminating data, such as cross-contamination from other samples in a multiplex sequencing run, can be identified by its low coverage. Scaffolds derived from the bacteriophage PhiX174 and Lambda DNAs that are frequently used as molecular standards in Illumina protocols can also be detected. Assemblies of yeast genomes with high heterozygosity, such as interspecies hybrids, often contain two types of scaffold: regions of the genome where the two alleles assembled into two separate scaffolds and each has a coverage level C, and regions where the two alleles co-assembled (collapsed) into a single scaffold that has a coverage level 2C Visualizing the data with Coverage-vs.-Length (CVL) plots, which can be done using Microsoft Excel or Google Sheets, provides a simple method to understand the structure of a genome assembly and detect aberrant scaffolds or contigs. We provide a Python script that allows assemblies to be filtered to remove contaminants identified in CVL plots.

The reported point centromeres of Scheffersomyces stipitis are retrotransposon long terminal repeats.

Point centromeres, found in some ascomycete yeasts such Saccharomyces cerevisiae, are very different in structure from the centromeres of other eukaryotes. They are tiny and nonrepetitive and contain only two short conserved sequence motifs. Until recently, point centromeres were thought to have a single evolutionary origin, in the budding yeast family Saccharomycetaceae. Most yeasts outside this family have centromeres that are many kilobases in size. Some have centromeres consisting of a large inverted repeat sequence, others have centromeric clusters of retrotransposons, and a third group including Candida albicans has centromeres with no conserved sequence features. It was recently reported that Scheffersomyces stipitis has point centromeres with a strongly conserved 125-bp core sequence, which is unexpected because S. stipitis is only distantly related to the known point-centromere species. We show here that the 125-bp core sequence is actually part of the long terminal repeat (LTR) of the Ty5-like retrotransposon Tps5, which forms a cluster in the centromeric region of each S. stipitis chromosome. Thus, the LTR of a centromere-associated retrotransposon confers centromere-like mitotic stability when cloned into a plasmid. The centromeric regions of S. stipitis contain three types of Tps5 element (Tps5a, Tps5b, and Tps5c) and a noncoding nonautonomous large retrotransposon derivative.

Population genomics shows no distinction between pathogenic Candida krusei and environmental Pichia kudriavzevii: One species, four names.

We investigated genomic diversity of a yeast species that is both an opportunistic pathogen and an important industrial yeast. Under the name Candida krusei, it is responsible for about 2% of yeast infections caused by Candida species in humans. Bloodstream infections with C. krusei are problematic because most isolates are fluconazole-resistant. Under the names Pichia kudriavzevii, Issatchenkia orientalis and Candida glycerinogenes, the same yeast, including genetically modified strains, is used for industrial-scale production of glycerol and succinate. It is also used to make some fermented foods. Here, we sequenced the type strains of C. krusei (CBS573T) and P. kudriavzevii (CBS5147T), as well as 30 other clinical and environmental isolates. Our results show conclusively that they are the same species, with collinear genomes 99.6% identical in DNA sequence. Phylogenetic analysis of SNPs does not segregate clinical and environmental isolates into separate clades, suggesting that C. krusei infections are frequently acquired from the environment. Reduced resistance of strains to fluconazole correlates with the presence of one gene instead of two at the ABC11-ABC1 tandem locus. Most isolates are diploid, but one-quarter are triploid. Loss of heterozygosity is common, including at the mating-type locus. Our PacBio/Illumina assembly of the 10.8 Mb CBS573T genome is resolved into 5 complete chromosomes, and was annotated using RNAseq support. Each of the 5 centromeres is a 35 kb gene desert containing a large inverted repeat. This species is a member of the genus Pichia and family Pichiaceae (the methylotrophic yeasts clade), and so is only distantly related to other pathogenic Candida species.

TPP riboswitch-dependent regulation of an ancient thiamin transporter in Candida.

Riboswitches are non-coding RNA molecules that regulate gene expression by binding to specific ligands. They are primarily found in bacteria. However, one riboswitch type, the thiamin pyrophosphate (TPP) riboswitch, has also been described in some plants, marine protists and fungi. We find that riboswitches are widespread in the budding yeasts (Saccharomycotina), and they are most common in homologs of DUR31, originally described as a spermidine transporter. We show that DUR31 (an ortholog of N. crassa gene NCU01977) encodes a thiamin transporter in Candida species. Using an RFP/riboswitch expression system, we show that the functional elements of the riboswitch are contained within the native intron of DUR31 from Candida parapsilosis, and that the riboswitch regulates splicing in a thiamin-dependent manner when RFP is constitutively expressed. The DUR31 gene has been lost from Saccharomyces, and may have been displaced by an alternative thiamin transporter. TPP riboswitches are also present in other putative transporters in yeasts and filamentous fungi. However, they are rare in thiamin biosynthesis genes THI4 and THI5 in the Saccharomycotina, and have been lost from all genes in the sequenced species in the family Saccharomycetaceae, including S. cerevisiae.

Evolutionary instability of CUG-Leu in the genetic code of budding yeasts.

The genetic code used in nuclear genes is almost universal, but here we report that it changed three times in parallel during the evolution of budding yeasts. All three changes were reassignments of the codon CUG, which is translated as serine (in 2 yeast clades), alanine (1 clade), or the 'universal' leucine (2 clades). The newly discovered Ser2 clade is in the final stages of a genetic code transition. Most species in this clade have genes for both a novel tRNASer(CAG) and an ancestral tRNALeu(CAG) to read CUG, but only tRNASer(CAG) is used in standard growth conditions. The coexistence of these alloacceptor tRNA genes indicates that the genetic code transition occurred via an ambiguous translation phase. We propose that the three parallel reassignments of CUG were not driven by natural selection in favor of their effects on the proteome, but by selection to eliminate the ancestral tRNALeu(CAG).

Ploidy Variation in Kluyveromyces marxianus Separates Dairy and Non-dairy Isolates.

Kluyveromyces marxianus is traditionally associated with fermented dairy products, but can also be isolated from diverse non-dairy environments. Because of thermotolerance, rapid growth and other traits, many different strains are being developed for food and industrial applications but there is, as yet, little understanding of the genetic diversity or population genetics of this species. K. marxianus shows a high level of phenotypic variation but the only phenotype that has been clearly linked to a genetic polymorphism is lactose utilisation, which is controlled by variation in the LAC12 gene. The genomes of several strains have been sequenced in recent years and, in this study, we sequenced a further nine strains from different origins. Analysis of the Single Nucleotide Polymorphisms (SNPs) in 14 strains was carried out to examine genome structure and genetic diversity. SNP diversity in K. marxianus is relatively high, with up to 3% DNA sequence divergence between alleles. It was found that the isolates include haploid, diploid, and triploid strains, as shown by both SNP analysis and flow cytometry. Diploids and triploids contain long genomic tracts showing loss of heterozygosity (LOH). All six isolates from dairy environments were diploid or triploid, whereas 6 out 7 isolates from non-dairy environment were haploid. This also correlated with the presence of functional LAC12 alleles only in dairy haplotypes. The diploids were hybrids between a non-dairy and a dairy haplotype, whereas triploids included three copies of a dairy haplotype.

Comparative genomics of biotechnologically important yeasts.

Ascomycete yeasts are metabolically diverse, with great potential for biotechnology. Here, we report the comparative genome analysis of 29 taxonomically and biotechnologically important yeasts, including 16 newly sequenced. We identify a genetic code change, CUG-Ala, in Pachysolen tannophilus in the clade sister to the known CUG-Ser clade. Our well-resolved yeast phylogeny shows that some traits, such as methylotrophy, are restricted to single clades, whereas others, such as l-rhamnose utilization, have patchy phylogenetic distributions. Gene clusters, with variable organization and distribution, encode many pathways of interest. Genomics can predict some biochemical traits precisely, but the genomic basis of others, such as xylose utilization, remains unresolved. Our data also provide insight into early evolution of ascomycetes. We document the loss of H3K9me2/3 heterochromatin, the origin of ascomycete mating-type switching, and panascomycete synteny at the MAT locus. These data and analyses will facilitate the engineering of efficient biosynthetic and degradative pathways and gateways for genomic manipulation.

Centromeres of the Yeast Komagataella phaffii (Pichia pastoris) Have a Simple Inverted-Repeat Structure.

Centromere organization has evolved dramatically in one clade of fungi, the Saccharomycotina. These yeasts have lost the ability to make normal eukaryotic heterochromatin with histone H3K9 methylation, which is a major component of pericentromeric regions in other eukaryotes. Following this loss, several different types of centromere emerged, including two types of sequence-defined ("point") centromeres, and the epigenetically defined "small regional" centromeres of Candida albicans Here we report that centromeres of the methylotrophic yeast Komagataella phaffii (formerly called Pichia pastoris) are structurally defined. Each of its four centromeres consists of a 2-kb inverted repeat (IR) flanking a 1-kb central core (mid) region. The four centromeres are unrelated in sequence. CenH3 (Cse4) binds strongly to the cores, with a decreasing gradient along the IRs. This mode of organization resembles Schizosaccharomyces pombe centromeres but is much more compact and lacks the extensive flanking heterochromatic otr repeats. Different isolates of K. phaffii show polymorphism for the orientation of the mid regions, due to recombination in the IRs. CEN4 is located within a 138-kb region that changes orientation during mating-type switching, but switching does not induce recombination of centromeric IRs. Our results demonstrate that evolutionary transitions in centromere organization have occurred in multiple yeast clades.